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Abstract 

We provide analytical and numerical evidence of the existence of classically stable, 
string-like configurations in a 2+1 dimensional analog of the Skyrme model. The 
model contains a conserved topological charge usually called the baryon number. 

Our strings are non-topological solitons which have a constant baryon number per 
unit length. The energy per length containing one baryon is, however, less than the 
energy of an isolated baryon (radially symmetric “baby Skyrmion”) in a region of 
the parameter space, which suggests a degree of stability for our configurations. In a 
limiting case, our configuration saturates a Bogomolnyi-type bound and is degenerate 
in energy per baryon with the baby Skyrmion. In another limiting case, the energies 
are still degenerate but do not saturate the corresponding Bogomolnyi-type bound. 
Nonetheless, we expect the string to be stable here. Both limiting cases are solvable 
analytically. 

The Skyrme modeP is a non-linear sigma model containing topological solitons which de¬ 
scribe the low-energy dynamics of mesons and baryons. Its 2-1-1 dimensional analog was 
studied in various contexts. Here the fields take values in the two-sphere 5^, hence the 
dynamical variables correspond to maps from to S^. Imposing that all configurations 
go to a constant at spatial infinity effectively compactifies the spatial also into an 5^, 
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and hence all maps are characterized by a topological charge corresponding to the winding 
number of maps from 3“^ to With the addition of the Hopf term to the action, it was 
shown that the solitons have fractional spin and statistics^. The stability and dynamics 
of the solitons of the model, called baby Skyrmions, have been studied^ and it was shown 
that stable localised solitons exist with the incorporation of a Skyrme term and a mass 
term. For a study of many solutions of the CP^ model in 2+1 dimensions, see reference [4]. 
Soliton-soliton scattering has been studied in reference [5] which gives a simpler analog of 
the corresponding phenomenon in the usual Skyrme model®. In a condensed matter context 
it was shown that 2+1 dimension Skyrmions are also relevant in quantum Hall systems and 
anyonic superconductors^. 

The energy functional is given by the usual kinetic term of the 0(3) non-linear sigma 
model®, a four derivative term which is analogous to the Skyrme term®’®, and a mass term 
(or interaction with an external magnetic held) which actually breaks the symmetry and 
picks out a vacuum®®, 

E = ^ J cPx di${x) ■ di(j){x) + {di(j){x) x d 2 (p{x))‘^ + /i^(n — 0(T))^ . (1) 

Here a = 1,2,3 are the components of a unit vector 0, ^(T) ■ 0(T) = 1, and n is a 
constant unit vector taken for convenience to be (0,0,1). The kinetic term along with the 
Skyrme term are not sufficient to stabilize a baby Skyrmion contrary to the usual Skyrme 
model. The kinetic term in 2+1 dimensions enjoys (suffers from) conformal invariance and 
the baby Skyrmion can always reduce its energy by inhating (inhnitely). Hence one adds 
the mass term which limits the size of the baby Skyrmion. The usual Skyrme term of course 
prohibits the collapse of the soliton. 

The conhguration giving rise to a baby Skyrmion with topological charge N is (see the 
hrst article of reference [2]) 

(j){r,6) = (sin/(r) cos iV6*, sin/(r) sin iV6*, cos/(r)) (2) 

where /(r) goes from tt at the origin to zero at +cx), and where the topological charge, which 
by analogy, we call the baryon number, is given by 
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( 3 ) 


B = ^ [ (fxB{x) = ^ f (f'x (f ■ di(j){x) X d 2 (p{x). 

4:71 J 471 J 

This gives (for TV = 1) the ordinary, non-linear differential equation for /(r)^: 


r -|- 


sin^ f{r) \ fit 


/"(r)+ 1 


sin2/(r) I /'(r)sin/(r)cos/(r) 
-.2 X - 


( 4 ) 


\f\r) 

\ J \ J 

sin f(r) COS f(r) 2 • £ / \ r\ 

-— ixj. _ sin j- i^r) = 0. 

This equation can only be integrated numerically, except in some limiting cases. A solution 
with unit baryon number B, has /(O) = tt, descending with a finite slope at the origin 
and arriving to zero (mod 27i) at infinity with an exponential, two-dimentional Yukawa- 
like fall-off governed by the mass term. The solution is singular in the sense that it is not 
differentiable at the origin, however, the energy density is well defined everywhere. 

To construct the string-like configuration, we simply re-interpret the radial coordinate 
r and the angular coordinate 6 as two Cartesian coordinates: {r,d) {x,y). This has the 

effect of laying the configurations which previously occured along rays, from the origin to 
r = -|-cx), in a linear progression along the y axis with the value (0,0, —1) at a; = 0 and 
the value (0, 0,1) (the vacuum) at a; = -|-cx). This way we obtain one baryon per length 
271 in the y direction for the right half plane. Along the y axis the configuration has the 
value (0,0,—1) which is not in the vacuum direction. We must extend our configuration 
into the left half plane in a way that we also reach the vacuum at x = —oo. The symmetric 
fashion of achieving this is to glue on to the line at a; = 0 a configuration which corresponds 
to creating the same string-like configuration as before, however, after having performed a 
rotation by tt. This rotation has the effect of reversing the directions of x and y. Indeed the 
configuration, where we have rescaled y, 


1( \ i ■ f( \ 7r|/sign(a;) . . 7iysign{x) 

(p{x, y) = (sm J [X) cos---, sm / (a;) sm---, cos j (a;)) 


( 5 ) 


L ^ ' I 

with /(O) = TT, /(±cx)) = 0 mod 27i, has unit baryon number per length L and has a baryon 
number density which is independent of y. 

7rsign(a;) . 


B{x,y) = 


smf{x)f'{x). 


( 6 ) 


3 



The function f{x) satisfies the differential equation 

(l + (0 sin^ f{x)^ f{x) + (^0 sm2f{x) _ ^2 ^ q_ 

The solution is found in the right half plane on imposing the boundary condition /(O) = tt, 
f{+oo) = 0. The continuation to the left half plane is done by either reflecting /(—|a;|) = 
/(|a;|) so that f{—oo) = 0 or continuing smoothly the solution which then interpolates from 
/(O) = TT to f{—oo) = 271. Either interpolation gives the same energy and baryon number 
density. Equation (7) can actually be integrated analytically giving rise to the quadrature: 



Equation (8) is not terribly useful to obtain f{x), since we still must invert the function 
dehned in terms of the integral on the left hand side, however, we may use it to obtain the 
energy per baryon as the following integral: 

rl ( 47^2 2 \i /2 

= 8T dy{l + ^(1 - yVj ^ ' (9) 

This allows for a numerical calculation of the energy per baryon without recourse to numer¬ 
ical resolution of any differential equation. Subsequently we minimize with respect to L to 
hnd the actual minimum energy string conhguration. We hnd Estring = 1.55356 x An with 
Lstring = 3.4542 for the value of the mass parameter /i^ = 0.1 as chosen in the second article 
of reference [3] whereas for the baby Skyrmion Eskyrmion = 1.564 x 47r. Comparing Estring 
to Eskyrmion (obtained numerically using the shooting method) for general we hnd the 
following curve for Estring — Eskyrmion (see Figure 1). This means that the string will be 
stable against disintegration into individual baryons in the region 0 < /i < 1. The upper 
limit, /i = 1, is obtained numerically with an error of less than 0.002%, but we lack an 
analytical understanding of this fact. The string with baryon number N is obtained from 
(5) by the straightforward substitution y Ny. It is evident from the structure of the 
ansatz that a string with baryon number N has length NLgtring and energy NEgtring- 
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Taking the limit /i —0 while rescaling the coordinates appropriately, removes the Skyrme 
and mass terms, and we obtain the conformally invariant 0(3) non-linear sigma model with 
energy 

E = ^ J (fx di(j){x) ■ di(j){x). (10) 

As is well known, this can be written in a form making the Bogomolnyi bound^^ apparent 


E = - J X ((9i0(a;) ± 0 x d 2 (j){x)] ± 2B{x) . 


Hence E > 47r|-B|, the Bogomolnyi bound is saturated when 


(9i0(a;) ± 0(a;) x d 2 (t){x) = 0, 


depending on the sign of the baryon number. This equation has a solution for the string 
conhguration, the equation of motion being 

= 0 (13) 

with solution 


f{x) = 2 arctan^CK 


where a is an arbitrary scale parameter. The solution (14) shows that each half string 
becomes inhnitely wide achieving tt ai x = -|-cx) for the -|- sign, and at a; = —cx) for the — 
sign. Equation (12) also has a baby Skyrmion type solution, satisfying the equation 


f’{r) = ± 


sin /(r) 


and given by 


/(r) = 2arctan (ar)^^ . 


The N baryon generalization also satishes the Bogomolnyi bound. Since the string saturates 
the Bogomolnyi bound per baryon, it is degenerate per baryon with a conhguration of an 
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isolated N baryon solution that also saturates this bound. Hence the string conhguration is 
classically stable in this limit. To elaborate this further we simply calculate the energy for a 
conhguration 0 = 00 + <50, where 0o satishes the equation (13) or (15), using the expression 


( 11 ) 

E(0) = E(0o) + /d2f 
= E(0o) + SE. 


01(50 ± (00 X 02^0 + (50 X 0200 + (50 X 02(50) 


( 17 ) 


The energy of the huctuation 6E is clearly a positive, semi-dehnite quantity. 

The other limit, essentially /r —+cx) also yields an interesting and analytically solvable 
system. Here, if we scale the coordinates appropriately while taking the limit, we can 
dispense with the kinetic term, leaving the energy fnnctional 


E 


1 

2 


cPx 


{Slip X d2$f + y(n 



(18) 


The eqnations of motion with the ansatz (2) and (5) are eqnally well integrable. Indeed for 
the radially symmetric sitnation we get 


yt y y Zi 


(19) 


which integrates to 


cos/(r) = —- ij;r‘ + 1, 


( 20 ) 


with the bonndary conditions /(O) = tt, and /(r) —0 as r becomes large. The interesting 
featnre of this solntion is that the held achieves (0, 0,1) exactly at a hnite radius i?, ][E) = 0. 
We hx R by minimizing the energy which amonnts to the condition that the conhguration 
achieves (0, 0,1) smoothly and yields R = 2/yjji. We glue on to the ontside for r > i? the 
vacnum. The energy turns ont to be 
4 

Eskyrmion g dvTp.. (21) 

For the N baryon ansatz (2)), this energy is simply multiplied by N. For the ansatz (5), we 
obtain the following equation of motion 
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which also is trivially integrated, 


— sin f{x) = 0 


(22) 




71 I 2 


71 / 


(23) 


exhibiting the same, intriguing behaviour that the soliton cuts off at a hnite half-width 
W = 27r/(/iL). The energy per baryon is also given by 
4 


Estring r, ^71 fl. 


(24) 


Interestingly enough, this does not depend on L, the length per baryon, and the energy 
per baryon of a thin string is equal to that of a wide string. The area per baryon is also 
independent of the value of L 


A 


string 


4:77 

= 2 X L xW = —. 

/i 


(25) 


Amazingly enough this corresponds exactly to the area occupied by the radially symmetric 
baby Skyrmion, Askyrmion = = 47r//i. We have not been able to unearth the reason 

for this apparent equality of areas, and also for the apparent shape invariance of the string 
itself. There is perhaps an underlying area preserving diffeomorphism invariance in the 
model resulting in this “incompressibility”. 

This limiting case does in fact contain a Bogomolnyi type bound, which is not saturated. 
Indeed, 


E = - j (fx ((9i0 X (920 ± /i(n — (9))^ ± 2jx^B{x) 


(26) 


and the extra cross term in the perfect square involving n integrates to zero since it is a 
total divergence. Hence the energy per baryon satishes 


E > 471 


(27) 


In any case, since our strings per baryon number are degenerate with an isolated baryon, 
and we actually do not expect a lower energy configuration in this sector, our configurations 
seem to be stable against disintegration into single baryons. 
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The classical stability analysis for general /x is not straightforward. For fi G [0,1] we 
have found that the string is energetically stable against disintegration into single baryons. 
This is not valid for disintegration into conhgurations with B = 2. Taking the spherical 
ansatz for this “baby deuteron”, n = 2 in (2) and solving (numerically) the corresponding 
differential equation for /(r), yields a conhguration that is quite tightly bound, 

Edeuteron gj^y^Ynion 

For example with = 0.1 we get A = 0.191 x 47r, as found in the second article of reference 
[3]. Hence, if not classically, we expect a string to be able to quantum mechanically tunnel 
into a conhguration of isolated lumps with B = 2. Calculation of the tunelling rate requires a 
detailed understanding of the potential on the space of conhgurations, which for the moment 
is lacking. In the two limiting cases studied above, the baby deuteron is degenerate with 
the string (taken in lengths with baryon number two) and we expect that it will be stable 
even quantum mechanically (although in the second case there may be some other bound 
conhguration for the deuteron). 

As there is no topological reason for the stability of our strings, even less so than for 
usual cosmic strings which do enjoy and utilise the topology in the transverse direction for 
part of their stability, an inhnite string will reduce its energy per unit length (to zero) by 
diluting itself (inhnitely) along its length. This is not the relevant consideration. We expect 
our conhgurations to be important when considering the low energy excitations of lumps 
with large, but hnite baryon number, N. The topology of indicates that the minimum 
energy conhguration will be a localised, lump of energy, in a held conhguration with most 
likely some large symmetry group. The excitation of this system will always preserve the 
baryon number. Hence if a string conhguration is excited, the relevant energy consideration 
must be done at hxed baryon number. This boils down to considering the energy per baryon, 
as we have done above. Excitation of this matter will involve the string conhgurations, be 
they only meta-stable. This should correspond to a diherent low energy phase for the matter 
described by this energy functional. 
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